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ABSTRACT 

We report the discovery of 33 Lyman break galaxy (LBG) candidates at z 8 detected in Hubble 
Space Telescope (HST) Wide Field Camera 3 (WFC3) imaging as part of the Brightest of Reionizing 
Galaxies (BoRG) pure-parallel survey. The ongoing BoRG survey currently has the largest area (274 
arcmin^) with Yogs (or Fios), -^125, and i^ieo band coverage needed to search for z ~ 8 galaxies, about 
three times the current CANDELS area, and slightly larger than what will be the final CANDELS 
wide component with Yi^^ data (required to select z ~ 8 sources). Our sample of 33 relatively 
bright lo98-dropout galaxies have J125 band magnitudes between 25.5 and 27.4 mag. This is the 
largest sample of bright ( J125 ^ 27.4) z ~ 8 galaxy candidates presented to date. Combining our 
dataset with the Hubble Ultra-Deep Field (HUDF09) dataset, we constrain the rest-frame ultraviolet 
galaxy luminosity function at z ~ 8 over the widest dynamic range currently available. The combined 
datasets are well fitted by a Schechter function, i.e. ^(L) = (j)^{L/L^)^ without evidence for 

an excess of sources at the bright end. At 68% confidence, for h = 0.7 we derive (j)^ = (4.3I21) ^ 
Mpc-^, = -20.26^^-3^, and a very steep faint-end slope a = -l.9Stoil- While the best-fit 
parameters still have a strong degeneracy, especially between (j)^ and M^, our improved coverage at 
the bright end has reduced the uncertainty of the faint-end power-law slope at z ~ 8 compared to the 
best previous determination at ±0.4. With a future expansion of the BoRG survey, combined with 
planned ultradeep WFC3/IR observations, it will be possible to further reduce this uncertainty and 
clearly demonstrate the steepening of the faint-end slope compared to measurements at lower redshift, 
thereby confirming the key role played by small galaxies in the reionization of the universe. 
Subject headings: cosmology: observations — galaxies: evolution — galaxies: formation — galaxies: 
high-redshift 



1. INTRODUCTION 

Finding the earliest galaxies in the universe and char- 
acterizing their properties and contribution to the reion- 
ization of the universe are some of the most import goals 
of extragalactic astronomy. The Wide Field Camera 3 
(WFC3) aboard the Hubble Space Telescope (HST) has 
significantly expanded the high-redshift frontier with the 
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detection and study of galaxies at z > 7. The ultradeep 
WFC3/IR observations of the Hubble Ultra-Deep Field 
(HIJDFOQ) and nearby fields have so far yielded 73 z ^ 7 
and 59 z ^ 8 Lyman- break gal axy (LBG) candidates 



f )Ouwens et al. 2Qllb| see also 
cLure et al.||201l| ), including o 



Lorenzoni et al. 

nr 



one at z 



sno 



2011 



et al.||2Ullap. THese ultradeep observations show tan- 
talizing evidence for a rapid evolution of the rest-frame 
ultraviolet (UV) galaxy luminosity function (LF) from 
z = 6 to z = S and a declining star-format ion rate with 
increasing redshift ( Bouwens et al.|2Qlla|b ) , as expected 
on the basis of dark matter halo assembly (Trenti et al. 
|2QlQp . 

However, because the ultradeep datasets cover only a 
few WFC3/IR fields, their dynamic range is limited by 
the small volume that they probe. In particular they pro- 
vide poor constraints on the population of bright g alaxies 
at z ^ 8, wh ich are very rare and highly clustered ( Trenti 
et al.||2Q12b|). A complete understanding of the number 
density and overall shape of the LF at z ~ 8 requires a 
large-area survey to search for the brightest galaxies at 
these epochs. 

Identifying the brightest z ^ S candidates from broad- 
band photometry is also of fundamental importance to 
provide the best targets for spectroscopic follow-up stud- 
ies aimed both at confirming the redshift of the sources 
as well as inferring the properties of the intergalac- 
tic medium (IGM) in proximity of the se galaxies (e.g., 
Schenker et al. 2012 Treu et al. 2012). Spectroscopic 
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confirmation of tiie redsiiift of LBG galaxies has been 
carried out for lar^e samples at z ~ 4 — 6 (e.g., Malhotra 



et al.|2QQ5|| Stark et al.f2Qll|) , 
to z = T.^l pno et al.p OliD. 



and recently extended out 

_ These studies show that 

photometrically selected samples of LBGs have very low 
contamination (~ 10%, see Malhotra et al. 2005). It is 
necessary to extend the frontier of spectroscopy further 
into the epoch of reionization, at z 8, not only to pro- 
vide definitive proof that the LBG selection continues to 
be reliable for candidates into the reionization epoch, but 
more importantly to infer the ionization state of the IGM 
from the study of the Lyg equivalent w idth distribution 
of LBG sources (e.g., |Treu et al.||2012 ). 

At present there are several ongoing WFC3/IR surveys 
designed to identify a relatively large sample of bright 
(^AB ^ 27; Mab ^ —20) z ~ 8 galaxies, namel y the 
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CANDELS Mu l ti-Cycle Treasury (MCT) program (jGro- 
gin et al.||2011l iKoekemoer et al.pOll), the Hubble in- 



frared Pure Parallel Ima ging Extragalactic Survey (HIP- 
PIES) fYan et al. |2011b[), and the Brightest of Reionizing 
Galaxies (BoRG) ( ['Irenti et al.|2QTH |2012b[ ). In the near 
term, the sample size of bright z ^1 galaxies will also 
be augmented by cluster le nsing surveys (e.g.. Hall et ah] 
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HUDF09, 



2012| [Bradley et al.||2012[ ), including the Cluster Lens- 
ing And Supern ova survey with Hubble (CLASH; Post- 



man et al.|2012l) MCT program, which recently reported 
the discovery of a magnified LBG candidate at z ~ 9.6 
(|Zheng et al.|2012|). All of these programs take advantage 
of mult i- band HST optical NIR data to search for z ~ 8 
galaxies as Y^c^s, or ^105 dropo uts using the Lyman-break 



technique ( Steidel et al.|p!996 ) 

Currently, the shape ot the bright-end of the galaxy 
UV LF at z = 8 is debated. Theoretical and nu- 
merical investigation s predict that it should remain 
Schechter-like (Schechter 1976), of the form (j){ L) = 
(t ■ 
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Figure 1. Comparison of area and depth of current surveys that 
can search for z ~ 8 LBGs. Because only ~ 260 arcmin^ of the 
planned CANDELS wide survey includes F-band imaging (red and 
orange), the ongoing BoRG survey (blue; ~ 274 arcmin^) provides 
the largest area for a y-band dropout survey. The depth of our 
BoRG survey is only ~ 0.3 mag shallower than the wide part of 
the CANDELS survey. The ultradeep HUDF09 WFC3/IR obser- 
vations are illustrated by the green point. 



sample of lo98 dropouts described in 



Trent i et al. 



(2011 



at z - 7 - 10 (Trenti et al. 



2010 



2012b) to include an additional 139 arcmin^ of area. Our 
latest catalog now contains eight very bright (> 8cr) 
z ~ 8 candidates and an additional 25 bright 2: ~ 8 
candidates detected at lower significance (> 5cr), with 



Jaacks et al.||2012D. Conversely, some numerical studies ^ ^ ^z=8' This is currently the largest-area search for 



have suggested a possible excess of sources at the bright 
end, with a non-Schechter behavior tied to the inefficient 
onset of AGN feedback at early times (e.g., Finlator et al. 
201 Ip . " 

Qbservationally, two recent papers from the CAN- 
DELS dataset in GOOD S-South yielded different re- 
sults. Yan et al. (2011a) found a significant excess of 
8 candidates with Mab < —21.0 over the 



bright z 

expectation from a Schechter function using data from 
- 80 arcmin^ of area in the GOODS-South field. The 
shape of their z ^ 8 luminosity function resembles a 
step function (see their Figure 4). On the other hand, 
analyzed the same dataset and de- 
Schechter function for bright z ^ 8 



Oesch et aL| (|2012| 
rive a well-behave( 
sources, consistent with the extrapolation of the LF evo 
lution fr om z ^ 7 and with the pr edictions of the LF 
model of [Trenti et"aL] ( |2010| ). The [Oesch et al^ ( |2012 ) 
result (M*^8 = -19.8f_o[lj) is also consistent with the 
first-epoch BoRG (BoRG09) determination of the knee 
of the Schechter function (M*_8 = -20.2 ± 0.3, with the 
error bar at a fixed a and (j)^). Our initial BoRG data 
were about 0.5 mag shallower than the current CAN- 
DELS dataset, but we already had larger sky coverage 
(~ 140 arcmin^ vs. ~ 80 arcmin^). 

Here we take advantage of the additional fields recently 
acquired as part of the BoRG survey and update the 



z ^ ^ candidates, now totaling ~ 274 arcmin^. Addi- 
tionally, given the random pointing nature of our pure- 
parallel HST program, the BoRG dataset is disti nct in 
that it is m inimally affected by cosmic variance dTrentil 
I& Stiavelli^ 2008). Therefore, this catalog is uniquely po- 
sitioned to set the tightest constraints on the number 
density of the brightest z ^ 8 galaxies and the bright 
end of the rest-frame ultraviolet galaxy luminosity func- 
tion {Muv < -19.6; L > O.SLJ^g)- 

This paper is organized as follows. We describe the 
BoRG survey in Section [2] We present the observations 
and photometry in Section [3] and dropout selection in 
Section [U The results and constraints on the z ^ 8 UV 
LF are discussed in Section (H We summarize our results 
and conclusions in Section |6] Throughout this work, we 
assume a cosmology with = 0.3, = 0.7, and Hq = 
70 km s~^ Mpc~^. This provides an angular scale of 4.8 
(proper) kpc arcsec"^ at 2: = 8.0. All magnitudes are 
expressed in the AB photometric system ( Oke|[1974 ). 



2. THE BORG SURVEY FOR z - 8 GALAXIES 

In 2009 we initiated the BoRG survey, a pure- 
parallel WFC3 survey that complements deep and ultra- 
deep WFC3/IR observations by looking for very bright 
(rriAB < 27; Mar < -20) galaxies at z > 7.5 ([Trenti 



et al. [2011 , 2012b) by obtaining WFC3 imaging in tour 
filters (1^^606 W, F098M, F125W, F160W) on random 
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Table 1 

B0RGO9 Survey Fields, Exposure Times, and 5cr Limiting Magnitudes^ 





Field 


Alternate 


FID 


Q:J2000 


<^J2000 


F600LF 


F606W 


F098M 


F125W 


F160W 






Name^ 




(deg) 


(deg) 


t (s) 




t (s) 


7T7, J A ™ 


t (s) 


7T7, J A ™ 


t (s) 


7T7, ? A ™ 


t Cs) 


7T7, 7 A ™ 


borg 


^.0214+1255 


yan24 


11702 


33.410 


12.915 


2294 


25.9 






2806 


25.9 


1403 


25.9 


1403 


25.7 


borg 


^_0540-6409 


borg2n 


11700 


84.879 


-64.151 






3171 


26.9 


4112 


26.6 


2309 


26.6 


1406 


26.3 


borg 


^-0553-6405 


borg81 


11700 


88.276 


-64.088 






3624 


27.0 


6418 


27.0 


2612 


26.8 


2012 


26.3 


borg 


5_0624-6432 


borg2t^ 


11700 


95.898 


-64.528 






2133 


26.7 


1806 


26.3 


1206 


26.4 


503 


25.6 


borg 


^.0624-6440 


borg2k 


11700 


95.951 


-64.663 






2135 


26.7 


2909 


26.6 


1206 


26.6 


906 


26.0 


borg 


^_0637-7518^ 


borg93 


11700 


99.265 


-75.313 






4290 


26.7 


6218 


26.7 


2412 


26.6 


1612 


26.0 


borg 


5_0751+2917 


borgOt 


11700^ 


117.709 


29.282 


3732 


26.6 


2826 


26.8 


18641 


27.4 


5115 


27.1 


3912 


26.8 


borg 


^_0756+3043 


borgOc 


11700 


118.989 


30.718 






2600 


26.7 


4712 


26.7 


1906 


26.6 


1406 


26.2 


borg 


^_0808+3946 


borgln 


11700 


122.089 


39.759 






2600 


26.5 


4612 


26.4 


2206 


26.6 


1406 


26.0 


borg 


5_0819+4911 


borgOg 


11700 


124.830 


49.184 






1908 


26.4 


3009 


26.6 


1206 


26.5 


806 


25.8 


borg 


^.0820+2332 


borg30^ 


11700 


125.014 


23.535 






2556 


26.7 


3109 


26.5 


703 


26.1 


703 


25.8 


borg 


5_0906+0255 


borg73 


11700 


136.405 


2.925 






3106 


26.9 


5518 


27.0 


2709 


27.0 


1906 


26.6 


borg 


5_0909+0002 


borg66 


11700 


137.286 


-0.030 






2650 


26.8 


3909 


26.8 


1806 


26.7 


1006 


26.0 


borg 


^.0914+2822 


borg39 


11700 


138.569 


28.362 






2571 


26.8 


4615 


26.8 


2206 


26.8 


1706 


26.5 


borg 


^_0922+4505 


borglr 


11700 


140.406 


45.088 






2708 


26.6 


4812 


26.6 


2106 


26.5 


1706 


26.3 


borg 


^.0926+4000 


borg45 


11700 


141.393 


40.006 






1276 


26.2 


2806 


26.3 


1106 


26.2 


903 


25.9 


borg 


^_0926+4426 


yan28 


11702 


141.382 


44.426 


2374 


26.6 






6012 


27.1 


1603 


26.7 


1403 


26.5 


borg 


^_1031+3804 


borg70 


11700 


157.715 


38.059 






1815 


26.4 


3109 


26.4 


1506 


26.3 


1306 


26.0 


borg 


5_1152+5441 


borgOy 


11700 


177.958 


54.684 






2898 


27.0 


6021 


27.0 


2809 


27.1 


1906 


26.7 


borg 


5_1153+0056 


borgOj 


11700 


178.182 


0.931 






2647 


26.8 


4515 


26.7 


2209 


26.7 


1606 


26.4 


borg 


^.1209+4543 


borgOp 


11700^ 


182.355 


45.724 


2234 


26.6 


2707 


27.0 


13729 


27.4 


3709 


27.2 


2909 


26.8 


borg 


^_1230+0750 


borglv 


11700 


187.470 


7.825 






2436 


26.6 


4112 


26.1 


1806 


26.0 


1406 


25.6 


borg 


5_1242+5716 


yanll 


11702 


190.554 


57.270 


2800 


26.5 






5215 


27.0 


2509 


26.9 


2309 


26.6 


borg 


^_1245+3356 


borg49 


11700 


191.186 


33.936 






1789 


26.7 


3409 


26.8 


1506 


26.7 


1106 


26.2 


borg 


^_1337+0028 


yanl9 


11702 


204.202 


-0.464 


2270 


26.4 






6818 


27.0 


1203 


26.6 


1203 


26.3 


borg 


5_1341+4123 


yan32 


11702 


205.131 


41.384 


3810 


27.0 






17435 


27.6 


3206 


27.2 


2806 


26.9 


borg 


5_1437+5043 


borg58 


11700 


219.234 


50.719 






2754 


26.9 


4912 


27.0 


2509 


27.0 


1806 


26.5 


borg 


^_1524+0954 


yanSl 


11702 


231.041 


9.906 


2078 


26.3 






8718 


27.0 


1603 


26.5 


1303 


26.2 


borg 


5_1632+3737 


borglk 


11700 


247.892 


37.609 






1260 


26.4 


2909 


26.7 


1206 


26.6 


906 


26.0 


Note. — The total survey 


area of BoRG09 is 


135 arcmin^, with an effective search 


area for lo98 dropouts 


of ~ 98 


arcmin^ 


after 



accounting for incompleteness. Unless otherwise noted, each field has an area of 4.7 arcmin^. 
^ 5cr magnitude limits i n a r = 0.32^^ aperture, corrected for Galactic extinction. 
^ Used in Trenti et aLl(|201l|). 

^ Field has multiple, partially overlapping exposures for a total area 6.9 arcmin^. Exposure times quoted are the sum of all exposures. 
^ Data missing due to scheduling constraint /conflict. 



^ Also includes data from HIPPIES program 11702. 

sighthnes at high Galactic latitudes {\b\ > 30°). Because 
luminous massiv e galaxies at these re dshifts are expected 
to be clustered ( Trenti et al.||2Q12b ), the random point- 
ing nature of our BoKG survey is ideal to mitigate the 
severe effects of c osmic variance ( Trenti fc Stiavelli|2QQ8 



place z ^ S galaxies within reach of ground-based spec- 
troscopic follow-up observations, as has been attempted 
on some of our brighter candidates ( Schenker et al.|2Q12 
Treu et al.f2Q12 l). 



[Robertson 2010). In fact, our current survey geometry 



of 59 independent WFC3/IR fields makes our number 
counts of z ^ S galaxies essentially follow a Poisson dis- 
tribution and can be used to constrain the luminosity 
function as well as a contiguous single-field survey with 
about two times more area at the same depth. If all 
the BoRG area had been in a single field of 16.5 x 16.5 
arcmin^, we derive from "Trenti & Stiavelli' (2008) that 
cosmic variance would have dominated over Poisson un- 
certainty for a sample of 33 candidates (22% vs. 17% 
respectively, for a total fractional error of 28%). 

By reducing the uncertainty in the number density 
of bright sources due to cosmic variance, we can place 
stronger observational constraints on (or M^), the 
characteristic luminosity of the LF. This, in turn, can 
help break the well-known degeneracy between and 
a, the faint-end slope, when fitting a Schechter (1976) 
luminosity function to the data. Placing tighter con- 
straints on the value of a at redshifts 2; > 6 is crucial in 
determining the contribution of galaxies to the reioniza- 
tion of the universe. Further, the relative brightness of 
the LBGs discovered in the BoRG survey can possibly 



Our ongoing BoRG survey has covered 274 arcmin to 
date thanks both to its continuation in Cycle 19 (GO 
12752, PI Trenti) and by assimilating data from the sim- 
ilar HIPPIES pure-parallel program ( |Yan et al.||20Tlbl) 
as well as the coordinated parallel observations acquired 
as part of the Cosmic Origins Spectrograph (COS) GTO 
program. As illustrated in Figure [l] the BoRG survey is 
~ 0.3 mag shallower than the "wide" part of the CAN- 
DELS program that has Yios-band data, but has a larger 
area with lo98 (or Y105), -^125, and i^ieo band coverage 
that is needed to search for z ~ 8 galaxies. 

In this paper we focus on z ~ 8 LBG candidates se- 
lected as lo98-band dropouts. To use a homogeneous 
selection of z ^ S candidates, we specifically exclude the 
HIPPIES Cycle 18 dataset, which replaced the F098M 
band with the F105W band. We leave the identification 
of z ^ 8 Yios-band dropouts for discussion in a future 
paper. 

3. OBSERVATIONS AND PHOTOMETRY 

The primary goal of the BoRG survey is to identify 
galaxies at z ~ 8 and measure the bright end of the 
LBG rest-frame ultraviolet LF at this redshift. Our 
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Table 2 

B0RGI2 Survey Fields, Exposure Times, and 5cr Limiting Magnitudes^ 





Field 


FID 


J zUUU 


J ZUUU 


F600LF 


F606W 


F098M 


F125W 


F160W 




















Z \b) 




Z \b) 




Z \b) 




borg 


^.0110-0224^ 


11700^ 


17.532 


-2.395 


3892 


26.7 


13911 


26.9 


39297 


27.0 


13538 


26.9 


10032 


26.6 


borg 


^_0228-4102 


11541 


36.987 


-41.026 


1414 


26.5 






3406 


26.9 


1403 


26.8 


1403 


26.5 


borg 


5_0240-1857 


11541 


40.114 


-18.954 


1360 


26.4 






3406 


26.8 


1403 


26.7 


1403 


26.5 


borg 


5_0427+2538 


11533 


66.690 


25.640 


1200 


24.1 






2409 


25.2 


703 


25.1 


503 


25.0 


borg 


^_0436-5259 


11520 


69.059 


-52.986 


1932 


26.6 






4406 


27.1 


2906 


27.3 


2206 


26.8 


borg 


5_0439-5317 


11520 


69.855 


-53.278 


1932 


26.6 






4206 


27.0 


3006 


27.2 


2306 


26.8 


borg 


5_0440-5244 


11520 


69.959 


-52.731 


1932 


26.5 






6609 


27.2 


2003 


27.0 


1403 


26.5 


borg 


^.0835+2456 


12025 


128.821 


24.936 






4698 


26.8 


4209 


26.7 


2206 


26.7 


2009 


26.3 


borg 


^_0846+7654 


11520 


131.593 


76.893 


1401 


26.4 






4406 


27.0 


2003 


27.1 


1603 


26.6 


borg 


5_1010+3001 


12025 


152.406 


30.018 






9142 


27.2 


5612 


26.8 


3212 


26.8 


2812 


26.6 


borg 


^_1014-0423 


11524 


153.523 


-4.379 


1221 


26.3 






1909 


25.6 


1106 


26.2 


703 


25.9 


borg 


^_1031+5052^ 


12025 


157.691 


50.862 






9991 


27.4 


9629 


26.9 


4812 


27.1 


4212 


26.8 


borg 


5_1033+5051^ 


12025 


158.212 


50.860 






9640 


27.3 


6412 


26.8 


3212 


26.8 


2812 


26.5 


borg 


5_1051+3359 


12024 


162.822 


33.985 






4220 


26.9 


6235 


26.8 


3318 


27.0 


1912 


26.4 


borg 


^_1103-2330 


12025 


165.808 


-23.506 






8558 


27.0 


9429 


27.1 


4412 


27.2 


3612 


26.8 


borg 


5_llll+5545 


12025 


167.737 


55.751 






6524 


27.1 


5518 


26.9 


2606 


27.0 


2409 


26.6 


borg 


5_1119+4026 


11519 


169.514 


40.398 


1308 


26.2 






2206 


26.6 


1403 


26.8 


1403 


26.4 


borg 


^_1131+3114 


11519 


172.876 


31.289 


1316 


26.5 






2106 


26.5 


1403 


26.8 


1403 


26.4 


borg 


^_1301+0000 


11702 


195.318 


-0.007 


2150 


26.1 






5812 


26.4 


1603 


26.1 


1303 


25.9 


borg 


5_1408+5503 


12572 


211.993 


55.056 






3324 


26.9 


5623 


26.9 


2612 


26.9 


2612 


26.6 


borg 


^_1510+1115 


12572 


227.537 


11.242 






5326 


26.9 


8423 


27.1 


3812 


27.1 


3812 


26.8 


borg 


^_1555+1108 


12025 


238.857 


11.132 






4753 


26.9 


5515 


27.0 


2909 


27.0 


2509 


26.7 


borg 


5_1632+3733 


11700 


248.074 


37.557 






2751 


26.8 


5115 


26.9 


2406 


27.0 


1806 


26.5 


borg 


5_1815-3244 


11533 


273.628 


-32.734 


1200 


22.3 






2509 


21.8 


703 


21.6 


403 


21.5 


borg 


^.2057-4412 


11530 


314.340 


-44.207 


2500 


26.4 






5009 


26.7 


1203 


26.4 


803 


25.9 


borg 


5_2132+1004 


11524 


323.062 


10.064 


1355 


26.2 






2409 


26.5 


1006 


26.4 


503 


25.7 


borg 


5_2155-4411 


11541 


328.812 


-44.177 


2130 


26.7 






5609 


27.0 


1403 


26.7 


903 


26.2 


borg 


^.2203+1851 


11534 


330.705 


18.850 


3200 


26.8 






17229 


27.4 


2006 


26.8 


2806 


26.8 


borg 


^_2345+0054 


11702 


356.261 


-0.902 


2028 


26.4 






5612 


27.0 


1403 


26.8 


1403 


26.5 


bors 


5_2351-4332 


11528 


357.650 


-43.525 


1050 


26.3 






11123 


27.4 


4209 


27.3 


2806 


26.7 



Note. — The total survey area of BoRG12 is ~ 139 arcmin^, with an effective search area for YoQS dropouts of ~ 115 arcmin^ 
after accounting for incompleteness. The combined survey area of BoRG09+BoRG12 is ~ 274 arcmin^, with an effective search 
area for Yogs dropouts of ~ 213 arcmin^. Unless otherwise noted, each field has an area of 4.7 arcmin^. 
^ 5cr magnitude limits in a r = 0.32^^ aperture, corrected for Galactic extinction. 

^ Field has multiple, partially overlapping exposures for a total area 14.8 arcmin^. Exposure times quoted are the sum of all 
exposures. 

^ Field has multiple, partially overlapping exposures for a total area 5.9 arcmin^. Exposure times quoted are the sum of all 
exposures. 

^ Also includes data from HIPPIES program 11702. 



Table 3 

Photometry of YoQS-dropout {z ~ 8) Candidates in the 8cr Catalog 



ID 



Ck:j2000 



^J2 



J125 J125- H160 S/Nv'' S/No98 S/N125 S/N160 



borg 


^.0440-5244.682 


69.9455843 


-52.7320162 


25.9 + 0.1 


> 2.1 


0.0 + 0.3 


-0.8 


-0.5 


9.1 


5.7 


borg 


^_0751+2917_229 


117.7141714 


29.2715323 


26.5 + 0.2 


> 2.5 


-0.0 + 0.3 


-0.9 


0.0 


8.9 


6.5 


borg 


^.0909+0002.595 


137.2731625 


-0.0297391 


26.1 + 0.2 


1.8 + 0.6 


-0.4 + 0.4 


-0.5 


1.8 


8.6 


3.2 


borg 


^_1033+5051_126 


158.1863098 


50.8416866 


26.0 + 0.2 


> 2.5 


-0.2 + 0.3 


0.0 


-0.6 


8.1 


5.6 


borg 


^_1301+0000_160 


195.3070838 


-0.0189297 


25.5 + 0.2 


2.1 + 0.6 


-0.4 + 0.3 


0.7 


1.7 


9.5 


5.2 


borg 


^.1437+5043.1137 


219.210672 


50.7260085 


26.1 + 0.1 


> 2.7 


0.0 + 0.2 


-1.5 


-1.0 


10.9 


7.9 


borg 


5_1555+1108_1417 


238.8651549 


11.1393576 


26.6 + 0.2 


1.9 + 0.7 


-0.6 + 0.4 


1.1 


1.4 


8.4 


3.4 


borg 


5.2203+1851.1061 


330.6930404 


18.8581986 


26.0 + 0.2 


2.5 + 0.6 


-0.3 + 0.3 


0.2 


1.6 


9.4 


6.6 



^ Total magnitudes (AUTOMAG). 

^ Vqoq whenever possible; otherwise VeooLP- 

observations are designed to acquire WFC3 imaging in 
four filters (F606W, F098M, F125W, F160W), which are 
obtained on random and discrete sighthnes inherent to 
this being a pure-parahel program. We identify bright 
(^125 ^ 27) high-redshift galaxies by searching for Iqqs 
dropouts using the well-kno wn Lyman-break dropout 
technique ( Steidel et al.|1996l ). The J125 and i^ieo bands 
are used for source detection and to measure the rest- 
frame UV color, while the WFC3/UVIS FeoG band is 
used to reject low-redshift {z ~ 1.5 — 2) interloper galax- 



ies, which are the primary source of contamination to our 
z > 7.5 sample. As discussed in Section [4j our BoRG 
survey is optimized to minimize the probability of con- 
tamination from both low-redshift interlopers and cool 
dwarf stars. 

We also included data fr om the similar H IPPIES 
WFC3 pure-parallel program Yan et al. ( |2Qllb ) and the 
coordinated parallel observations from the (JOS GTO 
program. The filter selection of the HIPPIES program 
and the earlier part of COS GTO program differs from 
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Table 4 

Photometry of logs-dropout (z ~ 8) Candidates in the 5cr Catalog 





iU 


Q:J2000 


()J2000 


7- a 
-'125 




^098— Jl25 


J125— ^160 


Q / AT b 




Q / AT-, 


C / /\r 
'^/^Vi60 


borg 


;_0436-5259_1233 


69.0303878 


-52.9717897 


27.1 ± 0. 


,2 


> 1.8 


-0.4 ± 0.5 


1.4 


0.9 


5.5 


2.6 


borg 


_0553-6405_4006 


88.2647181 


-64.0821631 


26.7 ± 0. 


,2 


> 2.3 


0.1 ± 0.3 


0.2 


0.6 


6.8 


4.7 


borg 


_0751+2917_920 


117.706444 


29.2977181 


27.0 ± 0. 


,2 


> 2.3 


-0.6 ± 0.4 


-0.3 


-0.0 


7.1 


3.3 


borg 


_0756+3043_437 


118.9794155 


30.7177854 


26.5 ± 0. 


.2 


> 2.2 


-0.0 ± 0.3 


1.3 


0.9 


6.8 


4.6 


borg 


_0835+2456_253 


128.8067291 


24.9267442 


26.2 ± 0. 


,2 


2.0 ± 0.8 


-0.2 ± 0.4 


-1.2 


1.1 


6.8 


4.1 


borg 


_1031+3804_213 


157.7102797 


38.0497128 


26.6 ± 0. 


,3 


> 1.9 


-0.1 ±0.4 


0.3 


0.7 


5.4 


4.0 


borg 


;_1031+3804_831 


157.7353387 


38.0673682 


26.6 ± 0. 


.3 


> 1.9 


-0.5 ± 0.5 


1.4 


0.3 


5.5 


2.9 


borg 


_1103-2330_1180 


165.7881655 


-23.4990798 


26.7 ± 0. 


,2 


> 2.2 


0.0 ± 0.3 


-0.3 


-0.4 


7.8 


5.8 


borg 


_1131+3114_1244 


172.8573532 


31.2942314 


26.2 ± 0. 


,2 


> 2.0 


-0.3 ± 0.3 


0.8 


-0.2 


7.7 


4.3 


borg 


_1152+5441_1087 


177.9750532 


54.6979452 


27.2 ± 0. 


.3 


> 1.8 


-0.2 ± 0.4 


1.2 


0.9 


5.5 


3.2 


borg 


_1153+0056_540 


178.1909931 


0.9320074 


27.0 ± 0. 


,3 


1.9 ± 0.8 


-0.5 ± 0.4 


0.2 


1.1 


5.9 


3.2 


borg 


_1242+5716_159 


190.5672023 


57.2567197 


26.4 ± 0. 


,2 


> 2.2 


-0.2 ± 0.4 


0.0 


0.4 


6.6 


4.3 


borg 


_1408+5503_749 


212.0126402 


55.0585147 


26.5 ±0. 


.2 


> 2.0 


-0.2 ±0.3 


-1.1 


0.8 


6.7 


4.4 


borg 


;_1408+5503_980 


212.0082405 


55.0672755 


27.0 ±0. 


,2 


1.9 ±0.8 


-0.4 ±0.4 


-0.2 


1.0 


6.0 


3.5 


borg 


;_1437+5043_172 


219.2223469 


50.7080907 


27.1 ±0. 


,2 


2.0 ±0.8 


-0.4 ±0.5 


1.0 


1.0 


5.8 


2.8 


borg 


;_1437+5043_879 


219.2240496 


50.7259683 


27.3 ±0. 


.3 


> 1.8 


-0.3 ±0.5 


0.3 


0.5 


5.0 


2.7 


borg 


;_1510+1115_51 


227.5348783 


11.2225448 


26.7 ±0. 


,2 


1.9 ±0.8 


-0.2 ±0.4 


1.2 


1.0 


6.2 


3.9 


borg 


;_1510+1115_1236 


227.5521577 


11.2522441 


27.2 ±0. 


,3 


> 2.0 


-0.4 ±0.4 


1.1 


0.7 


6.3 


3.5 


borg 


;_1510+1115_1404 


227.5425951 


11.2615405 


26.7 ±0. 


,2 


> 1.9 


-0.3 ±0.4 


0.8 


0.7 


6.5 


4.0 


borg 


;_1555+1108_595 


238.8429322 


11.1279017 


27.3 ±0. 


,3 


> 1.9 


-0.2 ±0.4 


0.8 


-0.1 


5.7 


3.2 


borg 


;_1555+1108_1166 


238.8438689 


11.1420588 


27.2 ±0. 


,2 


> 1.9 


-0.3 ±0.4 


1.2 


0.2 


5.8 


2.9 


borg 


;_1632+3733_694 


248.0628393 


37.5568592 


27.4 ±0. 


,3 


> 2.0 


-0.4 ±0.5 


0.9 


0.7 


6.2 


2.7 


borg 


;_2132+1004_24 


323.0575947 


10.0443562 


26.5 ±0. 


,2 


> 2.1 


0.0 ±0.4 


1.4 


0.6 


5.3 


3.2 


borg 


;_2155-4411_341 


328.8301017 


-44.1819169 


26.6 ±0. 


.2 


> 2.5 


-0.4 ±0.4 


0.0 


-1.1 


7.3 


3.4 


borg 


;_2155-4411_1192 


328.8031162 


-44.1737506 


26.9 ±0. 


,3 


> 2.1 


-0.1 ±0.5 


-0.7 


0.6 


5.2 


3.1 



Note. — Sources reported in the 8a catalog (Table [Sf are not duplicated here, but of course are included in the 5cr catalog. 
^ Total magnitudes (AUTOMAG). 
^ Vqqq whenever possible; otherwise Vsqolp- 
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Table 5 

Candidates in the BoRG58 Protocluster ( [Trenti et al.|2012b| > 



ID 


<^J2000 


<^J2000 


J125" 


^098— -^125 


J125— ^160 


S/N6O6 


S/No98 


S/N125 


S/N160 


borg_1437±5043_1137 


219.210672 


50.7260085 


26.1 ±0.1 


> 2.7 


0.0 ±0.2 


-1.5 


-1.0 


10.9 


7.9 


borg_1437±5043_879 


219.2240496 


50.7259683 


27.3 ±0.3 


> 1.8 


-0.3 ±0.5 


0.3 


0.5 


5.0 


2.7 


borg_1437±5043_757 


219.2310489 


50.7240585 


27.1 ±0.2 


> 1.8 


-0.7 ±0.6 


1.1 


-0.5 


5.0 


1.7 


borg_1437±5043_435 


219.2202746 


50.7156344 


27.4 ±0.3 


> 1.8 


0.0 ±0.4 


0.6 


0.6 


4.9 


3.3 


borg_1437±5043_172 


219.2223469 


50.7080907 


27.1 ±0.2 


2.0 ±0.8 


-0.4 ±0.5 


1.0 


1.0 


5.8 


2.8 



Total magnitudes (AUTOMAG). 



Table 6 

Photometry of Other Possible Ybgs-dropout (z ~ 8) Candidates 



ID Q;j2000 (^J2000 -^125^ ^098" -^125 Jl25— ^160 S/Ny^ S/Nqqs S/N125 S/Niqq 

borg_0240-1857_392^ 40.0998766 -18.9604896 26.2 ± 0.2 2.1 ± 0.7 -0.1 ± 0.3 0.6 1.4 8.1 6.3 

borg_1632±3737_386d 247.8986483 37.6047539 25.1 ±0.1 1.9 ± 0.4 0.0 ± 0.2 -0.0 3.1 14.2 9.3 

^ Total magnitudes (AUTOMAG). 

^ Vqoq whenever possible; otherwise VeooLP- 

^ Unresolved source; possible QSO candidate or most likely an L/T dwarf star with an unusually red Yogs — -^125 color. 

^ Photometry contaminated by an adjacent bright source with similar colors {Yqqs — J125 = 1-6), but clearly detected in Veoe and Yqqs] 

most likely a low-redshift contaminant. 



BoRG in that those datasets used the F600LP filter in- 
stead of F606W to control for contamination fr om red- 
dened low-redshift sources. 'Trenti et al." ("2011) discuss 
the benefits of F606W compared to F600LP. Here we re- 
iterate that, because F606W has a larger transmission 
efficiency integrated in frequency space over the pass- 
band, it reaches deeper at fixed integration time. This 
also becomes apparent in the data, for example in field 
borg_0751 + 2917 (Table [T]), which has deeper data in 
F606W {miim = 26.8 with t = 2826 s) compared to 



F600LP {miim = 26.6 with t = 3732 s) despite the 
shorter integration time. 

Pure-parallel observations come with a unique set of 
challenges. Because the primary observations are spec- 
troscopic in nature, dithering is mostly absent in our 
pure-parallel WFC3 data. As a result, detector artifacts 
(e.g., hot/warm pixels), uncorrected cosmic rays, and 
WFC3/IR detector persistence must be carefully consid- 
ered in the data reduction and data analysis. To mitigate 
detector hot/warm pixels and cosmic rays, we employed 
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F098M 



F125W 



F160W 




1 borg_ 0909+0002_595^ 

F606W 



borgj 033+505 1_1 26 



Figure 2. Postage-stamp cutout images of four of the very bright (> 8cr) high-redshift Ybgs-dropout candidate galaxies. The cutout 
images are A" x A" ^ corresponding to 19.2 kpc on a side at z = 8, and are shown with a P.A.= 0°. 



a robust algorith m based on a v ariation of Laplacian 
edge detection (van Dokkum 2001 ) to filter the individual 
FLT files prior to comb ining them with MultiDrizzle 
(Koekemoer et al. 2002). The parameters for the algo- 
rithm were chosen to remove most sharp-edged artifacts 
from the data, but to also be rather conservative so as 
not to remove a ny real sourc es. 



As detailed in |Trenti et al. (2011), our observing strat- 
egy is also unique because it is designed to minimize the 
impact of WFC3/IR detector persistence from observa- 
tions taken in orbits immediately prior to our exposures. 
We take advantage of the property that the amount of 
residual image persistence appears to decay roughly as a 
power law with time. Briefly, observations in either the 
Ji25 band or i^ieo band are preceded in the same orbit 
by a comparably long lo98 band exposure, and when- 
ever possible, Ji25-band exposures precede those in i^ieo 
banq^ Therefore any persistent residual images from 
prior observations will be brighter in the bluer bands 
and thus will not contaminate our high-redshift dropout 

sample . 

MultiDrizzle (jKoekemoer et al. 2002[ ) was used to 
combine individual exposures in a given filter to pro- 
duce the final science images as well as the associated 
inverse- variance weight maps. The images are drizzled 



This Phase II strategy is implemented in GO 11700 and GO 
12752 (PI: Trenti), as well as the more recent COS-GTO parallel 
observations. 



to a final pixel scale of 0.08'Ypixel. Because BoRG is 
a pure-parallel program, the image depths vary among 
our fields. For a typical 4-orbit parallel field, we obtain 
exposures of 2200 s, 3800 s, 1800 s, and 1800 s in Veoe, 
^098, -^125, and i^i6o, respectively. In a {)A" diameter 
aperture, the corresponding 5cr limiting magnitudes are 
approximately 27.4, 27.1, 26.9, and 26.6, respectively. In 
Table [ij we list the survey fields, exposure times, and 5cr 
limiting magnitu des (r = 0.32^^ ape rture) of the initial 
BoRG09 dataset ( |Trenti et alT]|2011[ ), and in Table [2] we 
list the same properties lor the new fields in the signifi- 
cantly expanded BoRG12 dataset. We have released the 
reduced drizzled science and RMS images for the BoRG 
dataset, which is publicly available via the BoRG web- 
site at wolf359.colorado.edu and the High Level Science 

Products on the MAST websit^ 

We used SExtractor ( [Bertin fc Arnouts| [ l996[ ) in 
dual-image mode for object detection and photometry. 
For the detection image we used the WFC3/IR J125 
data. The inverse- variance weight images produced by 
MultiDrizzle were used to generate RMS maps for each 
image. We subsequently normalized the RMS maps to 
account for cor related noise introdu ced by the drizzling 
procedure (see Trenti et al.||2011l) and used them in 
SExtractor for both source detection and photometry. 
Sources were required to have at least 9 contiguous pix- 
els, each detected at a threshold of 0.7<j above the back- 
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F098M 



F125W 



F160W 




F600LP 



borg_1301+0000_160 



F606W 




borg_2203+1851_106 



Figure 3. Postage-stamp cutout images of four of the very bright (> 8cr) high-redshift 
images are A" x A" ^ corresponding to 19.2 kpc on a side at z = 8, and are shown with a P. A 



: 0°. 



-dropout candidate galaxies. The cutout 



ground. Object colors and signal-to-noise values were 
measured using isophotal magnitudes (ISOMAG) and 
total magnitudes were measured within scalable Kron 
apertures (AUTOMAG). 

The photometry in each fiel d has been correcte d for 
Galactic extinction using the Schlegel et al. (1998) dust 
maps. All measurements given in this paper include this 
correction, which is typically modest Ay < 0.2 because 
the survey primarily targets lines of sight at high galactic 
latitudes. 

4. SELECTION OF z ~ 8 yo98-BAND DROPOUT 
CANDIDATES 

We searched for 2: > 7.5 galaxies using a lo98-dropout 
selection criteria in two broadband colors. The general 
criteria are a strong break between the Iqqs and J125 
filters, a relatively blue or fiat J125 — i^ieo color, and 
a non-detection in t he optical Vma {o"^ VrnviLp ) band. 



Specifically, following [Trenti et al.| ( |2Qll[ |2Q12b[ ) we re- 
quire: 



S/Nvhand<^'^ 

(>098 - ^125) > 1.75 

( J125 - Hiqq) < 0.02 + 0.15[(ro98 - J125) - 1 



.75]. 



These criteria select galaxies with redshifts in the range 
z ~ 7.4 — 8.8 (see Section[5|. To minimize the probability 



of contamination of our sample by low-redshift interlop- 
ers, we impose a conservative non-detection threshold 
of 1.5(7 on the optical-band data. The non-detection in 
the optical V band is fundame ntal to produce a clean 
sample of z ~ 8 candidates (see Bouwens et al.|[2Qllb ). 
If a source is not detected in a given tilter, we set its 
magnitude to the corresponding la upper limit to calcu- 
late its colors. For our final catalog, we require sources 
to be detected with a signal-to-noise (S/N) threshold of 
6'/A^ > 5 in J125 and S/N > 2.5 in i^ieo, as measured in 
isophotal apertures (ISOMAG). For the bright 8a cata- 
log, we require sources to be detected with a signal-to- 
noise threshold of S/N > 8 in J125 and S/N > 3.0 in 

^160- 

Using these color criteria, we identified 33 relatively 
bright YoQS-dropout galaxy candidates with observed 
J125 band magnitudes between 25.5 and 27.4 mag over 
our 274 arcmin^ search area. Eight of these z ^ 8 LBG 
candidates are brighter than 26.6 mag and are detected 
at a significance level of > 8cr in J125. This is the largest 
sample of bright ( J125 ^ 27.4) z ^ 8 galaxies presented 
to date. We list the properties of the candidates in the 
8(7 and 5(7 catalogs in Tables [3] and [4j respectively. The 
postage-stamp cutout images of the eight 8a candidates 
are shown in Figures [2] and [3j 

The (lo98 - -^125) and ( J125 - i^ieo) colors of our z ^ 8 
candidates along with the color-color selection criteria 
are illustrated in Figure [4j We also show the expected 
colors of galaxies simulated over a wide range of redshifts 
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-0.8 -0.6 -0.4 -0.2 0.0 0,2 0.4 0.6 0.£ 

(•^125~^16o)^B 



Figure 4. lo98 — -^125 vs. J125 — i^ieo two-color diagram used 
to select our Ybgs-band dropout candidates. The colors of our 8cr 
sources are shown by the black data points. The error bars and 
lower limits are la (68% confidence). The gray region represents 
the lo98 — -^125 and ji25 — ^160 colors of our selection criteria. The 
colored points represent the expected colors of galaxies simulated 
over a wide range in redshifts. The g reen region in dicates the 
colors of low -mass L,T dwarf stars (e.g.,|Knapp et aI72004 Ryan 
|et al.||2bll ' Holwerda et al. 2012 (in prepjj . in addition to the 
color-color selection shown here, Yogs -dropouts also need to be 
undetected (S/N < 1.5) in the deep optical imaging performed in 
the V band by the BoRG survey. The non-detection in V band is 
fundamental to produce a clean sample of z ~ 8 candidates (see 
[Bouwens et al. 2011b). 



and also those of l ow-mass dwarf stars (e.g., |Knapp et al. 
[2004 ; Ryan e t al.||2011t Holwerda et al. 2012 (in prepjj. 
As l^^igure [4] shows, our conservative lo98 — -^125 > 1-75 
color criterion is effective in minimizing contamination 
from cool dwarf stars and low-redshift interlopers. 

4.1. Comparison with the Earlier BoRGOO Sample 

In this paper we employ an improved data reduc- 
tion that includes Laplacian filtering (see Section 3|). 
Therefore, we expect some photometric scatter within 
the measurement uncertainty with respect to the previ- 
ous ca talogs of lo98 dropouts published in [Trent i et al. 
( 2011[ 2012b). Indeed, this is the case: one of the four 
- ^ - - ^ "^(12011), source 
le catalog with 



bright ofopouts identified in Trenti et al. 
"BoRG Ik" , is now just marginally out of t : 
lo98 — ^125 = 1-70. This is not surprising, as that candi- 
date was at the edge of the selection window in the previ- 
ous photometry and we were already considering it likely 
{p ~ 60%) to be a contaminant that scattered i nto th e 
dropout selection (see Section 5.2 in Trenti et al.||2011|). 
While it was, and still is, uncertain whether this source 
is at z ~ 8, our artificial source recovery simulations (see 
Section [5| take photometric scatter into account when 
determining the effective volume of the survey, making 
the derivation of the luminosity function robust. 

With the new data reduction we also verified the pho- 
tometry for field BoRG58 (here borg_1437 + 5043), where 
we identified a z ~ 8 protocluster candidate [Trenti et al.| 
([2012bj). The improved photometry confirms the pre- 
vious measurements, although we note that two of the 
fainter z ^ 8 candidates have scattered out of our cur- 
rent catalog. In this new reduction and catalog, they are 




Redshift 



Redshift 



Figure 5. The S{z, m) magnitude-dependent redshift selection 
function of one representative BoRG survey field, borg_1437+5043, 
for both the 5cr (left) and 8cr (right) 2; ~ 8 source catalogs. These 
selection functions were obtained through simulations to recover 
artificial sources in the BoRG images, as discussed in Section fs] 
Through these simulations, we computed S{z,m) for each of the 
59 individual BoRG fields. 

detected in the J125 band at 4.98<j and 4.88<j, slightly be- 
low our formal threshold of 5cr detections (see Table [5|. 
It is worth noting that two new lo98-dropout sources 
are also detected in this field, just below the detection 
threshold, at S/N ~ 4.5 in J125, providing circumstantial 
evidence that the overdensity extends to fainter luminosi- 
ties in line with our theoretical and numerical predictions 
dTrenti et al.|2012b[ ). Deeper iJST imaging would be very 
useful to further investigate the nature of this overden- 
sity. 

4.2. Possible z ^ 8 Candidates Excluded from the Strict 

Sample 

While two additional candidates pass our 8a selection 
criteria (see Table [6|, we exclude them from our final 
catalog. One of these sources (borg_0240— 1857_392) is 
extremely compact and clearly unresolved. Given its 
compactness, the likely explanation is that source is 
a Galactic star with either a significant fluctuation in 
the photometric measurement or with an unusually red 
^098 — -^125 intrinsic color (at the ~ 2a level). Because 
this object is located right above the Galactic Center 
(/ = 354.36830034°, b = 23.48948421°), the source is 
most likely a reddened L or T dwarf star. An exciting, 
but much less likely, alternative would bea2:~8QSO, 
although the limited area of BoRG implies that this oc- 
curs with p < 5% based on the Willott et al.,(2010) QSO 
luminosity function predictions. 

The other source we excluded (borg_1632+3737_386) 
is extremely bright at 25.1 AB mag in the J125 band 
and located immediately adjacent (SExtractor extrac- 
tion flag "2") to an even brighter (J125 = 21.4 mag) fore- 
ground spiral galaxy. The spiral galaxy is tidally inter- 
acting with a second spiral and has colors similar to the 
dropout candidate (lo98 — -^125 = 1-2, ^125—^160 = 0.35, 
Veoe — J125 = 3.5), but it is clearly detected in Veoe- In 
pri nciple, the Y n98-dropout could have been lensed (e.g., 
see Wyithe et al. 2011), but because of the similar colors 
with the foreground source, we consider it much more 
likely that this dropout is part of the foreground sys- 
tem, with its Ji25-band photo metry partially co ntami- 
nated by emission-line flux (see Atek et al. 2011) and a 
^606 continuum that cannot be detected in the current 
data assuming its Vqqq — J125 is similar to that of the 
spiral galaxy. 

5. THE z^8 LBG LUMINOSITY FUNCTION 

We derive the completeness, C{m)^ and magnitude- 
dependent redshift selection function, S{z,m), of our 
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Figure 6. Effective comoving volume for our z ~ 8 selection as 
a function of J125 magnitude for both the 5cr (blue) and 8cr (red) 
source catalogs. The effective volume shown here takes into ac- 
count reductions from both photometric scatter and incomplete- 
ness. 
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Figure 7. The z ~ 8 galaxy UV LF from the current BoRG 
dataset. The stepwise LFs for our 5a and 8cr catalogs are shown 
in blue and green, respectively. The red points represent the step- 
wise LF derived by Bouwens et al. (2011b) for the ERS+HUDF09 
dataset. The blue line is the best-tit Schechter LF from combining 
the B0RG+ERS+HUDFO9 dataset, providing the widest dynamic 
range in luminosity that is currently available. Our results are con- 
sistent with a Schechter form of the UV LF and do not indicate an 
excess of bright 2; ~ 8 LBGs. 



dataset from simulations as described in lOesch et al.l 
(|2009l |2012|) . Briefly, artificial galaxies with a range of 
spectral energy distributions (SEDs), luminosities, red- 
shifts, and sizes are added to the real images. We then 
rerun our detection and selection procedure on the data 
in each individual field to determine both C{m) and 
S{z^m) (see Figure [5|. The simulations are based on 
z 4 galaxy images rescaled to the desired input magni- 
tude and to higher redshift using standard evolutionary 
relations. In particular we use a size scaling of (1 + ^)~^ 
as determined from LBGs in the range z ^3 — 7 (|Fer- 



uson et al. 2004 I Bouwens et al. 



201 op and adopt a 



I Bom 



20041 lUesch et a. 



ope distribution of 



contmuum si 

f:i = —2.5 ±0.4, motivated by the recent measurements 



Table 7 

BoRG Stepwise Determination of the 
z ~ 8 UV LF^ 

Muv (10-4 Mpc-3 mag-i)b 



-22.14 
-21.64 
-21.14 
-20.64 
-20.14 
-19.64 



< 0.015 
023+° °^^ 

U.UZO_Q Q-,^g 

n 1QS+0-106 
^•-L^o_0.072 

604+°-^^'^ 

1 296+°-^^^ 
4 504+4-366 



Note. — Assuming h = 0.7. 
^ For the 5cr sample. 
^ The errors are derived from the 68% 
Bayesian credible intervals for a Pois- 
son distribution. 



at z > 6 



Bouwens et al. 



Finkelstein et al. 



(2011); 



^ ( |2Q12 

Dunlop et al. (2012]). From tnese simulations we derive 
the ettective volume for our z 8 selection as a function 
of J125 magnitude as shown in Figure [6) By design these 
simulations take into account the reduction of the effec- 
tive area and volume probed by the survey because of 
photometric scatter and presence of foreground sources, 
including any persistence images in the Iqqs data. We 
also note that because we are selecting sources from the 
J125 images, the J125 — i^ieo color is not an unbiased 
measure of the rest-frame UV slope /3, so the simulations 
we perform are crucial to derive the effective volume of 
our selection, and show that we retain high efficiency, 
with S{z,m) > 0.8 (Figure |6|. 

We construct the stepwise luminosity function from our 
BoRG observations in 0.5 mag bins for the sample of all 
33 Fogs-dropouts and for the smaller, but more robust, 
sample of 8 sources with high S/N (J125 > 8cr) detec- 
tions. The results are shown in Table |7| and Figure [7| and 
take into account 34% contamination for the S/N > 8 
sample and 42% contamination for t he S/N > 5 sample 
based on improved estimates following [Trent i et al.|( |201l" , 
2012b). We take into account the contam ination levels 
of the legacy fields following O esch et al.| ([2012) and es- 
timate the BoRG contamination rate b y usm^ the Early 
Release Science (ERS) WFC3/IR data ( [Windhorst et al. 
|2011| ), which also used the F098M filter as the Y band, 
m combination with the GOODS F606W data degraded 
to match our relative Veoe and lo98 depths. The main 
limitation of this approach is the limited area of the ERS 
data, which is a factor 10 smaller than the BoRG area. 
Also, we assume that contamination does not vary as a 
function of luminosity for a given dataset, which is the 
working assumption of previous studies deriving the LF 
at high redshift, both in legacy and pure-parallel fields. 
In principle, the contamination might vary with source 
brightness, depending on the LF difference for contami- 
nants and 2; ~ 8 galaxies with the same near-IR colors. 
Without a spectroscopic follow-up survey, and given the 
lack of confirmed z ry 8 galaxies even in legacy surveys 
fSche nker et al.[[2012| [Treu et al.'2012), it is challenging, 
if at all possible, to mvestigate this in more detail. Yet, 
it is reassuring that the spectra of the two BoRG z ^ 8 
candidat es that have been obse rved at Keck with NIR- 
CAM by [Schenker et al.[ ([2012[) and [Treu et al.[ ([2012) 



allow us to exclude the scenario where these sources are 
z ^ 1.5 galaxies with strong emission lines, which are 
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Table 8 

Effect of the Contamination Fraction on the LF 
Determination 



1/(1-/) 


\nC 


(/)*(10-4 Mpc-3) 


M* 


OL 


1.0 


-24.56 


5.6 


-20.28 


-1.80 


1.25 


-23.96 


4.7 


-20.31 


-1.90 


1.5 


-23.84 


4.4 


-20.29 


-1.95 


1.73^ 


-23.98 


4.3 


-20.26 


-1.97 


2.0 


-24.35 


3.8 


-20.28 


-2.05 


2.25 


-24.80 


3.5 


-20.28 


-2.10 


2.5 


-25.31 


3.6 


-20.23 


-2.10 



Note. — Assuming h = 0.7. 
^ Our fiducial value for the contamination correction. 



Comparison of z ' 



Table 9 

' 8 LF Determinations in the Literature 



Reference 



log(/)* (Mpc-3) 



OL 



This Work 



Oesch et al. 


2012 


» 


Bouwens et c 


tn 


20 


lib) 


Liorenzom et al. (2 
Trent 1 et al. (2Ull 
'McLure et al. (2U1 


ol 


liouwens et £ 


ii. (i^UTDfc 



_o 07+0.26 
-0.29 
-17+0.40 

_o 90+0.74 
<J-^<J_0.27 

-3.0 
-3.4 (fixed) 

-3.46 
-2.96 (fixed) 



— 20 2fi"^^"^^ 

ZU.ZO_Q 34 

-20.10 ±0.52 
-19.5 

-20.2 ±0.3 
-20.04 (fixed) 

-19.5 ±0.3 



-1.98: 
-2.06^ 



2+0.23 
^-0.22 
-+0.45 
-0.37 

-1.91 ±0.32 
-1.7 (fixed) 
-2.0 (fixed) 
-1.71 (fixed) 
-1.74 (fixed) 



Note. — Assuming h = 0.7. 

exp ected to be the m ain contaminants of o ur selection 
(see |Atek et al ||2Q11| and [Trenti et alT]|2Qll ). 

To explore the ettect of contamination m BoRG, we 
performed the experiment of leaving the contamination 
fraction as a free parameter in the LF fit, as discussed 
below (see Table [s]). A different contamination frac- 
tion would, under our assumptions, simply scale up or 
down the BoRG LF measurement. Shifting up our data 
by ~ 0.2 dex, assuming no contamination, can be in- 
terpreted very conservatively as an upper limit to the 
bright end of the z 8 galaxy luminosity function. As 
Figure [t] shows, the low- and hi^-S/N samples provide 
a consistent determination of the luminosity function of 
lo98 dropouts in the BoRG dataset and the data in our 
faintest bins agree with (and are actually slightly lower 
than) the brightest bins in [Bouwens et al.| (2Qllb) from 
ERS and HUDF09 data. Finally, the measure of a strong 
clustering signal from the Bo RG survey argues a gainst a 
strong contamination fraction |Trenti et al. ( 2Q12b| ). 

We derive Schechter function parameters by perform- 
ing a maximum-likelihood fit to the data, assuming a 
Poisson distribution of the galaxy number counts in each 
magnitude bin. Specifically, we maximize the Poisson 
likelihood for observing N^^^ sources in a given magni- 
tude bin when N^^^ are expected to be observed from 
a given Schechter LF. The likelihood C is expressed as 
C = UjUi P(A^;j',A7,f ), where P{N^^^ , N^^'p) is the 
Poisson probability distribution and the products are 
taken over all the fields, j, and magnitude bins, i. To 
establish the best estimate of the overall shape of the UV 
LF, we combine our BoRG dataset with t hose f rom th e 
deeper ERS+HUDF09 data of Bouwens et al.| (|2011bl), 
providing the largest dynamic range in luminosity cur- 
rently possible. Given that the CANDELS observations 
yield very different preliminary results depending on the 



team t hat analyzed the data to search for z ~ 8 galax- 
ies (see Yan et ar]|2Qllal Oesch et al.||2Q12 ), we exclude 
those data from our analysis. 




Figure 8. la (68%) and 2a (95%) confidence intervals in the 
Schechter LF fit from BoRG±ERS±HUDF data (red fines) com- 
pared to the CANDELS±ERS±HUDF determination by Oesch] 
et al. (2012) (blue lines). The BoRG fit has a preference for a 
marginally brighter M* , but the two datasets are consistent within 
their 2a contours. In particular, we note that our M* versus a 
parameter values have much better constraints than the previous 
studies at z ~ 8. 



At 



(4.3l|^) X 



, for = 

-\ = - 



0.7 we derive 



20.26 



+0.29 
-0.34' 



and a very 



68% confidence 
10-"^ Mpc" 

steep faint-end slope a = — 1-981q"22- observed in Fig- 
ure [7| overall the best fit provides a very good description 
of the data. The covariance in the parameter values is 
smaller than that in the CANDELS-hERS-hH UDF anal- 
ysis by pesch et"aI] ( |2012D ( |Yan et al.| (|2011a) do not fit 
a LF to their data), but it is still signihcant, as shown 
in Figure m From this fi gure it is also clearly evident 
that our versus a parameter values have much better 
constraints than previous studies at 2; ^ 8 (Oesch et al. 
2012 ). Formally the bright end of the luminosity function 



has a marginally low value. However, this is compen- 
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Figure 9. Determination of the Schechter function faint-end slope 
a for the galaxy LF as a fun ction of redshift. Literature determina- 
ti ons are fr om Oesch et al.|( |2010] ) at 2; ~ 0.7 — 2.5 ( magenta), from 

fddy & Steidel (2UUy) at z ^ 'I — 3 f o range). from|Bouwens et al.| 
JU7 ) at z ~ 4 (cy an1, from[ Bouwens et al.| ( |2012| > at 2: ^ 5 - t) 
:een), and from Bouwens et al. (2Ullb) at 2; ~ '( (blue). Our new 
determination at z ~ 8 is shown in red. All error bars are Icr (68% 
confidence). Our latest determination provides some evidence that 
the LF is becoming steeper at 2; > 7, consistent with the 2; ~ 7 
measurement. 

sated by the low 0* such that the best-fit Schechter func- 
tion = X IQ-^ Mpc-^M^ = -20.1 ± 0.52) 



derived by [Bouwens et al" (2011b) is fully consistent with 
our deterniination at 68% confidence. As can be seen in 
Figure |7j the faintest three bins of the BoRG dataset 
match v ery closely the brightest bins in [Bouwens et al7 



([2011b"). The latest B oRG LF is al so very similar to our 
earlier fit from Trent i et al.[ (|2011[), where we had kept 
0* and OL fixed (see Table [9|. 

Table [S] shows the impact of varying the contamina- 
tion fraction / for the BoRG dataset and its associated 
contamination correction 1/(1 — /) that is applied to the 
intrinsic counts per magnitude bin to obtain the esti- 
mated observed counts. The table shows that the impact 
of contamination around our fiducial value / = 0.42 is 
relatively modest. We note that the highest likelihood 
is associated to a contamination fraction just below our 
fiducial value, providing indirect evidence to suggest that 
our derived contamination fraction is robust. Varying 
the contamination fraction around / = 0.33 has almost 
no impact on the determination of M*, as the LF is well 
sampled by the BoRG data in this luminosity range, but 
affects mostly a; if we have underestimated /, then the 
faint-end slope would be even steeper. 

As shown in Table [9j our determination of the z ~ 
8 luminosity function is consistent with previous work 
|McLure et al. [ [2010 



Lor enzoni et al . [[ 2 1 1 1 iTr ^iti et al 



2011 b| [UescFet al. 2012|, taking 



Oil I Bouwens et aJ 
mto account the significant uncertamties in all these mea- 
sures. Formally, our best fit prefers a low normalization 
for (\)^ and a brighter (similar to the z ~ 7 value), 
but the covariance between these parameters is very large 
(see Figure[9|. In particular our BoRG+ERS+HUDF LF 
fit provides a better constraint on the M*vs. a parame- 
ter uncertainties tha n the CANDELS-hE RS-hHUDF de- 
termination at z ~ 8 'Oesch et al.' (12012") . 
In comparing different datasets, it is important to note 



that the bright end is still sparsely sampled; for exam- 
ple, the brightest BoRG magnitude bin includes only one 
source. Therefore, small- number fluctuations are very 
significant and Poisson noise in the brightest bins can 
have a large impact on the fit. However, our maximum 
likelihood fit shows, independent of the versus (\)^ de- 
generacy, that the faint-end slope a is being constrained 
with growing accuracy. Of note, we reduce its la un- 
certainty at z ~ 8 to ±0.2 and show in Figure |9] that 
the steepening trend suggested by previous studies (e.g., 
[Bouwens et al . 2011b) is also supported by our LF fit. 

'I'he steep faint-end slope derived in our best-fit LF 
implies that faint galaxies are playing a key role as ma- 
jor producers of ionizing photons. For a Schechter LF, 
most of the luminosity density contri bution is at the fain t 
end (e.g., see left panel of Fig. 3 in |Trenti et al.[|2QTQ ). 
Furthermore, in the case of a ~ —2, there is a logarith- 
mically divergent contribution from sources below the 
detection limit of the survey. This implies that the to- 
tal ionizing flux produced by galaxies is sensitive to the 
exact measure of a (as well as on the extrapolation be- 
yond the detection limi t), as shown, for example, in the 
left panel of Fig. 4 of [Bouwens et aL (2012). Assum- 
ing that star formation continues to be emcient in lower 
luminosity galaxies down to Mab ^ —12 at z ^ 8 (the 



limit of atomic hydrogen cooli ng halos, e.g. see Trent i 



et al.|201Q Finlator et al.|20lT ), then we expect that the 



current tilST observations are only observing ~ 20% of 
the light present at z ~ 8. A direct proof that we are 
detecting only the tip of the iceberg of star formation 
during the epoch of reionization is provided by B.ST ob- 
servations of high-redshift, spectroscopically confirmed 
GRBs that have failed to detect host galaxies despite 
reaching u ltrafaint sensitivity {Mab ^ —17, see |Tanvir] 
et al.|2012[ ). The non-detection of galaxies at locations m 
the sky where it is confirmed that star formation is hap- 
pening (because a GRB explosion has occurred) implies 
that Ma b ^, — 15 galaxies w ere indeed the main ionizing 



sources CTrenti et al. 2012a), in full agreement with the 
interpretation of the Ll^' tit derived in this paper. 

Interestingly, our determination of the z ^ 8 luminos- 
ity function at the bright end is between the debated 
measurements from the CANDELS dataset obtained by 
( 2012| and |Yan et al., ( |2011a| ). 
a 



Oesch et al. 



However, 

our data clear ly ^ow a well-beh aved Schechter func- 
tion, similar to Oesch et al.| (2012]), w ithout the unusual 
shape derived by |Y an et al. ([2011 at). Our LF determi- 
nation and that by Oesch et al. ( |20l"2 ) are consistent at 
the ~ 2a level. The observed ditterences could be par- 
tially related to the env iro nment (i.e. cosmic variance 



fTrenti & Stiavelli'"2008)). [Oesch et al] ([2Q12| discuss 



that the ERS measurement by [bou wens et al.j ([2 011b 



could have been affected by an over density. However, 
given that the BoRG da taset is not affected by cosmic 

([2011b|, it 
iderde 



variance and agrees with [Bouwens et al 

well be that the CANDELIS-ISouth held is underdense, 



may 



which would explain the lower M^ value derived from 
that dataset. The upcoming BoRG observations and 
those in the CANDELS-North field wih help clarify the 
nature of the discrepancy and provide further improve- 
ments to the determination of the bright end of the z ^ 8 
luminosity function. 

6. SUMMARY AND CONCLUSIONS 
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We present the discovery of 33 Lyman break galaxy 
(LBG) candidates at z - 8 detected in HST /WFC3 
imaging as part of the Brightest of Reionizing Galax- 
ies (BoRG) pure-parallel survey. Our sample of bright 
lo98-dropout galaxy candidates have Ji25-band magni- 
tudes between 25.5 and 27.4 mag, obtained from 59 in- 
dependent lines of sight over a total area of 274 arcmin^. 
The pure-parallel nature of BoRG allows us to obtain 
an estimate of the galaxy luminosity function which is 
unaffected by large-scale structure uncertainty and dis- 
tinct from determinations using legacy surveys limited to 
a single, or a few, contiguous fields. 

With our new data we detect galaxies between — 22 < 
Mab ^ —19.75, demonstrating the bright end of the 
z ^ 8 luminosity function is well described by a Schechter 
form, similar to that found at lower redshifts. Our mea- 
surement of the number density of galaxies near the 
BoRG detection limit (Mab ^ T^^^ agreement 



with the results by Bouwens et al. (2011b). Their mea- 
surement is based on EKIS+HUDF09 observations that 
extend to much fainter luminosities, but they lack the 
area to detect the brighter galaxies at Mab ^ —20.5 
we identify in the BoRG survey. The two datasets are 
complementary and allow us to obtain the best determi- 
nation yet of the galaxy luminosity function at z ~ 8, 
with = (4.3t^j) X 10-"^ Mpc-^ = -20.261^;^^, 

and a = -1.98^0.22 fo^" ^ = O-^- 

Covariance in the parameters, albeit reduced thanks to 
the increase in the dynamic range of the fit, is still very 
significant. However, this affects primarily versus 
(^*. We find that our BoRG+ERS+HUDF LF fit pro- 
vides a better constraint on the versus a parameter 
uncertainties than the CANDELS+ER S+HUDF deter- 

The faint-end 



8 Oesch et al. (2012). 



mination at z 

slope uncertainty is starting to be reduced to a point 
where there is a hint of steepening compared to z < 6. 
This steepening is expected based on numerical and theo- 
retical models of galaxy formation in the epoch of reion- 
ization ( |Trenti et al.||2010| iJaacks et al.l|2012| ). As a 
consequence of this large abundance of faint galaxies, 
it is expected that such systems will dominate the to- 
tal star formation rate and ionizing photon production 



([Shull et al. 2012), which has recently been confirmed 
observationally by the non-detection of hos t galaxies in 
a sample of six z > b Ga mma- Ray Bursts (Trenti et al. 



"it^[20T2|. 

hter tit of th( 



Our determination of a steep 
e LF is robust against our es- 



2012a[ [Tanvir 
a for a Schechter 
timate of the contamination fraction of the BoRG sur 
vey (/ = 0.42 for the 5cr sample); as shown in Table Jsj 
if / were higher, a would be steeper. In the near fu- 
ture, the upcoming ultradeep observations of GO 12498 
(PI Ellis) will improve the determination of the number 
density of the faintest galaxies observable by HST, while 
scheduled BoRG and CANDELS-North observations will 
further increase the search area for the brightest galaxies 
at z ~ 8, allowing us to further tighten the constrains on 
the LF. 

The z ^ 8 candidates identified here are good targets 
for follow-up observations from ground and space ob- 
servatories. We have started a spectroscopic campaign 
to confirm their redshift, to measure the distribution of 
Lya equivalent width (which is related to the IGM ion- 
ization state), and also to rule out contamination from 



emission-line galaxies at z ^ 1.5 (iT reu et al.|[2QT2l ). The 
availability of multi-object spectrographs both on Keck 
(MOSFIRE) and Gemini South (FLAMINGOS-2), wih 
enable faster progress, especially in fields that have over- 
densities of sou rces similar to the protocluster candidate 
we reported in Trenti et al. ( [2012b} . In addition, the 
brightest BoRG sources are expected to have star forma- 
tion rates of 10 — 20 Mq yr~^ based on their rest-frame 
UV luminosity, making them prime targets for ALMA 
observations as they are expecte d to fall well within the 
telescope sensitivity based on the Carilli et al. (|2008|) pre- 
dictions. If these candidates are spectroscopically con- 
firmed, this would allow us to possibly extend the de- 
tection of [C ii] 158 /im, high-J CO emission lines, and 
perhaps dust emission as well, to z ~ 8, after the recent 
record established at z = 7.1 from a QSO host galaxy 
(Venemans et al.| |2012 ). Finally, Spitzer /IRAC observa- 
tions, combined with spectroscopic redshifts, would also 
be very useful to quantify the rest-frame optical proper- 
ties of these sources before the advent of the James Webb 
Space Telescope. 
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